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The three-flavor crystalline color-superconducting (CCS) phase of quantum chromodynamics 
(QCD) is a candidate phase for the ground state of cold matter at moderate densities above the 
density of the deconfinement phase transition. Apart from being a superfluid, the CCS phase has 
properties of a solid, such as a lattice structure and a shear modulus, and hence the ability to sus- 
tain multipolar deformations in gravitational equilibrium. We construct equilibrium configurations 
of hybrid stars composed of nuclear matter at low, and CCS quark matter at high, densities. Phase 
equilibrium between these phases is possible only for rather stiff equations of state of nuclear matter 
and large couplings in the effective Nambu-Jona-Lasinio Lagrangian describing the CCS state. We 
identify a new branch of stable CCS hybrid stars within a broad range of central densities which, 
depending on the details of the equations of state, either bifurcate from the nuclear sequence of stars 
when the central density exceeds that of the deconfinement phase transition or form a new family of 
configurations separated from the purely nuclear sequence by an instability region. The maximum 
masses of our non-rotating hybrid configurations are consistent with the presently available astro- 
nomical bounds. The sequences of hybrid configurations that rotate near the mass-shedding limit 
are found to be more compact and thus support substantially larger spins than their same mass 
nuclear counterparts. 



I. INTRODUCTION 

Matter in the interiors of neutron stars is compressed 
by gravity to densities which are by factors 5 to 10 larger 
than the density of an ordinary nucleus. At such densi- 
ties baryons are likely to lose their identity and dissolve 
into their constituents - deconfined quarks. It has been 
suggested long ago that quark matter may exist in the 
interiors of compact objects Within the deconfined 
phase unbound quarks form Fermi seas and attractive in- 
teractions bind them into Cooper pairs thus giving rise 
to quark superconductivity 0, y] ■ Along with the decon- 
finement phase transition the dense matter undergoes a 
second transition related to the restoration of chiral sym- 
metry (the dynamical symmetry of the strong interac- 
tion). This symmetry may remain broken in some color- 
superconducting phases, such as the color-flavor-locked 
phase. 

If compact (hybrid) stars featuring quark cores sur- 
rounded by a nuclear mantle exist in nature, they could 
provide a unique window on the properties of QCD at 
high baryon densities. Even if the relevant densities and 
sufficiently low temperatures could be reached in a lab- 
oratory in the future, the conditions prevailing in com- 
pact stars are different from those produced in acceler- 
ators; the matter is long-lived, charge neutral and in 
/^-equilibrium with respect to weak interactions. 

A central question of the theory of compact stars, 
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which we will address in this work, is whether the equa- 
tion of state of matter at high densities admits stable 
configurations of self-gravitating objects in General Rel- 
ativity featuring deconfined quark matter, and if so, are 
the gross parameters of these objects, like their mass and 
radius, compatible with the known astronomical bounds? 
The deconfinement phase transition from baryonic to 
quark matter leads to a softening of the equation of 
state which could lead to an instability towards a col- 
lapse into a black hole. The details of a stability analy- 
sis may depend on the theoretical model of deconfined, 
color-superconducting quarks. The Nambu-Jona-Lasinio 
(NJL) model is a non-perturbative low-energy approx- 
imation to QCD, which is anchored in the low-energy 
phenomenology of the hadronic spectrum. While dynam- 
ical symmetry breaking, by which quarks acquire mass, 
is incorporated in this model, it lacks confinement 
Previous studies of hybrid color-superconducting stars 
within the three-flavor NJL model displayed a general 
instability of hybrid stars towards collapse into a black 
hole [5]. Stable stars featuring two-flavor superconduct- 
ing (2SC) phases were obtained with typical maximum 
masses M ^ 1.7 Mq under a favorable choice of con- 
stituent quark masses (which are related to the value of 
the baryo-chemical potential at zero pressure) [fj] or by 
replacing the hard NJL cut-offs by soft form factors with 
parameters fitted to the same set of data 0, ll| • Heavier 
objects are obtained if a repulsive vector interaction is 
introduced in the NJL Lagrangian 

More phenomenological studies of hybrid stars based 
on the MIT bag model where carried out using a 
generic parameterization of the quark matter equation 
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of state [T3|- Non-perturbative QCD corrections to the 
equation of state of the Fermi gas were accounted through 
the parameter — I — c, where c ~ 0.3 is a reasonable 
estimate jTlj. For this value of the correction stable hy- 
brid stars containing color-flavor-lockcd superconducting 
quark matter exist with maximum mass 1.9 Mq. 

Because of (3 equilibrium in the light quark sector the 
chemical potentials of u and d quarks obey the constraint 
Hd — Hu + Me from which it is evident that the 

Fermi surfaces of the light flavor quark flavor are shifted 
by an amount corresponding to the chemical potential 
of electrons. Furthermore, the Fermi surface of strange 
quarks is mismatched with the Fermi surfaces of light u 
and d quarks because of the large strange quark mass. 
A generic feature of fermionic systems with mismatched 
Fermi surfaces is the existence of phases with broken spa- 
tial symmetries. A particular realization is the Larkin- 
Ovchinnikov-Fulde-Ferrell (LOFF) phase [li,[i3, which 
is theoretically predicted in different fermionic s yste ms 
ranging from dilute atom s |15| . nuclear matter (lal to 
dense quark matter [1, [13, [TMlia HI [HI [H, [H, [IJ 

The rotational and translation symmetries are broken 
by the superconducting LOFF phase, since the conden- 
sate carries a non-zero momentum ,13l The study of 
crystalline color superconductivity (CCS) in the context 
of QCD reveals that this phase could be the true ground 
state of deconfine d q uark matter at intermediate densi- 
ties d [13, H [S [la [m, [H . Later on we shall describe 
a particular realization of the three-flavor LOFF phase, 
which is used in our search for stable sequences of hybrid 
stars 0,13. 

The nuclear equation of state, as is well known, can 
be constructed starting from a number of different prin- 
ciples [H, We tested a large number of equa- 
tions of state to construct hybrid star configurations. 
These belong to the classes of (i) non-relativistic vari- 
ational and Bruckner-Hartree-Fock theories which use as 
an input a non-relativistic potential fitted to the elastic 
nucleon-nucleon scattering data; (ii) relativistic mean- 
field models which are fitted to the bulk properties of 
nuclear matter, and (iii) relativistic models which in- 
clude correlations at the level of the covariant scattering 
amplitude (Dirac-Bruckner-Hartree-Fock theories). For 
present purposes, it is sufficient to characterize these 
equations of states by their stiffness; as we shall see, only 
the stiffest equations of state are admissible for phase 
equilibrium between nuclear and quark matter. 

The paper is organized as follows. In Sec. Hi] we dis- 
cuss the input equations of state with the emphasis on 
the CCS phase of quark matter. Section IIIII contains 
our results for the sequences of non-rotating and rapidly 
rotating hybrid stars. Our conclusions are collected in 
SecHVl 



II. PHASE EQUILIBRIUM 
A. The three-flavor crystaUine phase of QCD 

The LOFF phase arises naturally in three-flavor QCD 
as a generic feature of a fermionic systems with mis- 
matched Fermi surfaces [l[il[iliii,|23,[2l|2i,[2l|23]. 
The simplest possible realization of the LOFF phase in- 
vokes a single plane-wave modulation of the order pa- 
rameter, related to a single value of the center-of-mass 
momentum of the Cooper pairs. More complex Ansatze 
for the order parameter maximize the free energy gain of 
the superconducting phase for periodic three-dimensional 
lattice structures (e. g. face-centered-cubic lattice in 
Ref. [1^ or body-centered-cubic lattice in Ref. [19]). The 
studies of complex structures can be carried out in the 
Ginzburg-Landau regime, where the order parameter is 
small, as is the case near the critical temperature of the 
phase transition or the critical mismatch for unpairing. 
Independent of the form of the lattice structure, the re- 
sulting CCS state breaks the translational and rotational 
symmetries since Cooper pairs always carry non-zero 
center-of-mass momentum with respect to some fixed ref- 
erence frame. Some of the characteristic features of the 
CCS state such as a non-zero shear modulus and gapless 
excitations imply that it could play a distinctive role in 
the phenomenology of compact stars [13, [iM ■ 

In the three-flavor case the quark condensate is written 
as a superposition of plane waves 

Pi 

(ipazC^bi'fij) (X^£ai3i<^iji^i^exp{2iqf' -r) . (1) 

1=1 m=l 

Here Greek (Latin) indices correspond to color (flavor); 
Pj is the number of plane waves of the particular crys- 
talline structure considered; Ai, A2, A3 refer to ds, us, 
ud pairing respectively, and 2q^ is the quark pair mo- 
mentum. The magnitude of q^„ and of the gap parame- 
ters A/ are determined by the minimization of the ther- 
modynamic potential fi. 

The studies of the crystalline structure of the LOFF 
phase revealed that there exists a window of the param- 
eter / fi, where Ms is the in-medium strange quark 
mass and /i the mean quark chemical potential, where 
two crystalline phases, called CubeX and 2Cube45z, 
could c omp ete for the ground state of three-flavor quark 
matter 'l9!|. For both of the structures Ai = 0, A2 = A3, 
P2 ^ Pz=P and |q2| = Iqaj. In the CubeX phase P = 4; 
for each pairing channel / = 2,3 the wave vectors {q} 
point to the vertices of a square, and the two squares are 
arranged in such a way that they point to the vertices 
of a cube. On the other hand, in the 2Cube25z one has 
P = 8; each wave vector set {q} forms a cube, and the 
two cubes are rotated by 45 degrees around an axis per- 
pendicular to one of the faces of the cube. The analysis of 
Ref. [1^ has been extended in Ref. [2^ by evaluating self- 
consistently the strange quark mass for each value of /z, 
therefore translating the window of into a window 



of /Lt. The LOFF phases then exist within the chemical 
potential range 442 MeV < ^ < 515 MeV j2fli]. Thus, 
having the dependence of the pressure on the chemical 
potential, we are in a position to construct phase equilib- 
rium and obtain the pressure as a function of the energy 
density. 

A straightforward normalization of the quark pressure 
in the NJL model requires that the pressure vanishes at 
zero density and temperature [2^[2^. In the terminology 
of the MIT bag model, this is equivalent to subtraction of 
a bag constant from the thermodynamic potential [loj . 
Since the value of the bag constant is related to con- 
finement which is absent in the NJL model, it appears 
reasonable to consider changes in its value, and hence in 
the normalization of the pressure. We shall consider the 
simple case of a constant shift in the asymptotic value of 
the pressure; alternatives include the use of form factors 
for the bag constant 01 or the Polyakov loop at nonzero 
density |,3ft. . 

In Ref. [13] the calculation of masses and condensates 
was carried out for rj — Gd/Gs = 0.75, where Gd 
and Gs are, respectively, the diquark and the quark- 
antiquark coupling constants, and the proportionality 
factor 1] is obtained via a Fierz rearrangement of the in- 
teraction term in the Lagrangian. This regime is usually 
referred to as "intermediate coupling" . Here we shall 
adopt a "strong coupling" regime with rj — 1 [l^, [sl] , for 
only in the latter case the matching to the nuclear equa- 
tions state can be performed without variations in the 
bag constant (we will discuss this point in more detail 
in the following subsection). Clearly, a more complete 
analysis will require a (re)computation of the effect of 
the strong coupling on the competition among the var- 
ious superconducting phases such as the CFL, 2SC and 
2-flavor LOFF phases; the possibility of multiple transi- 
tions between these phases is not considered here. 

For completeness we reproduce here the parameters of 
the NJL model of Ref. [20|] used in our study: the values 
of the bare m, d and s quark masses are m„ — rrid ~ 
5.5 MeV and = 135.4 MeV, the coupling constant 
for the 4-fermion term of the NJL Lagrangian is Gs = 
1.81/A^, the coupling constant for the t'Hooft term is 
K = 8.80/A5, where the ultraviolet cutoff is A = 643 
MeV. With these parameters the following observables 
of QCD are reproduced: the pion decay constant = 
93 MeV, the pion mass = 135 MeV, the kaon mass 
niK = 497 MeV, and the rj' mass m^- = 924 MeV. 



B. Matching the equations of state 

Physically, the true nuclear equation state must go over 
to some sort of quark equation of state at some density 
if deconfinement takes place in nature. Since we have 
only models of deconfined matter and nuclear matter, 
this transition is modeled by requiring that there exists 
a baryo-chemical potential at which the pressures of these 
phases are equal. This is equivalent to the condition that 
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FIG. 1: Pressure versus density for the models A {heavy, black 
online), Al {medium-light, red online), and B {light, blue on- 
line). The nuclear equations of state are shown by dashed 
lines {dark blue online). For the models A and Al the nuclear 
(low density) equation of state is the same; for the models A 
and B the quark (high density) equation of state is the same. 
At the deconfinement phase transition there is a jump in the 
density at constant pressure. 



pressure, P, vs. chemical potential, /j,, curves for these 
phases cross (matching). If the P(^) curves for chosen 
equations of state of nuclear and quark matter do not 
cross, the models are incompatible in the sense that they 
cannot describe the desired transition between nuclear 
and quark matter. The low-density equation of state of 
nuclear matter and the high-density equation of state of 
CCS matter are matched at an interface via the Maxwell 
construction. The phase with largest pressure is the one 
that is realized at a given chemical potential. Thus, ac- 
cording to the Maxwell construction of the deconfinement 
phase transition, there is a jump in the density at con- 
stant pressure as illustrated in Fig. [TJ The high-density 
regime contains two equations of states for crystalline 
color superconductivity which differ by the normalization 
of pressure at zero density (or, equivalently, the value of 
the bag constant). The model Al is normahzed such that 
the pressure vanishes at zero-density. For the models A 
and B the zero density pressure is shifted by an amount 
6p = 10 MeV/fm'^. As discussed above, this is equiv- 
alent to a variation of the bag constant whose value is 
related to confinement, is unspecified in the NJL model, 
and is uncertain in general. We are aware of the arbitrari- 
ness of the latter procedure, the sole practical purpose of 
which is to produce an equation of state which can be 
matched to a particular nuclear equation of state. Yet 
another possibility is to set Sp = 0, but vary the value 



of the constituent masses of the hght quarks in the fit 
of the parameters of the NJL model; for small values of 
the light quark masses the matching between quark and 
nuclear equations of state is facilitated Furthermore, 
we need to adopt a large value for the ratio 77, e. g. rj — 1, 
since for rj = 0.75 - the conventional value of this ratio 
- no matching occurs between the quark and the nuclear 
equations of state. 

A set of nuclear equations of states were tested for 
matching with the models above; it included about dozen 
equations of state, listed in Refs. jH,|23|. Only two equa- 
tions of state based on the Dirac-Bruckner-Hartree-Fock 
approach [2^ 1 are suitable to match with the quark equa- 
tions of state presented above. These are shown in Fig.[T] 
The selected equations of state are the two most hard 
equations of state in the collection. It should be noted 
that the matching with softer equations of state can be 
enforced by varying Sp by larger amounts than quoted 
above. For our purposes, however, the three equations 
of state that follow upon matching are sufficient. The 
values of the chemical potential (in MeV) and pressure 
(in MeV/fm'^) at the quark- nuclear matter interface are 
Ifi = 1234.2; p = 96) for the model Al, (1230.0; 95) for 
the model A and (1234.8; 108) for the model B. 



III. RESULTS 

A. Non-rotating configurations 

In this subsection we consider equilibrium and stability 
of cold hybrid stars with CCS cores; rotating configura- 
tions are discussed in the next subsection. Each equation 
of state defines a sequence of equilibrium, non-rotating 
stellar configurations in General Relativity, which can 
be parameterized in terms of the central density pc of 
the configuration. It is assumed that the configurations 
are cold, so that the temperature (or entropy) does not 
play any significant role in determining the stellar struc- 
ture. The spherically symmetric solutions of Einstein's 
equations for self-gravitating fluids are given by the well- 
known Tolman-Opennheimer-Volkoff equations (33j . A 
generic feature of these solutions is the existence of a 
maximum mass for any equation of state; as the cen- 
tral density is increased beyond the value corresponding 
to the maximum mass, the stars become unstable to- 
wards collapse to a black hole. One criterion for the sta- 
bility of a sequence of configurations is the requirement 
that the derivative dM/dpc should be positive (the mass 
should be an increasing function of the central density) . 
At the point of instability the fundamental (pulsation) 
modes become unstable. If stability is regained at higher 
central densities, the modes by which the stars become 
unstable towards the eventual collapse belong to higher- 
order harmonics. For configurations constructed from a 
purely nuclear equation of state the stable sequence ex- 
tends up to a maximum mass of the order 2 Mq (Fig. [2]); 
the value of the maximum mass is large, since our chosen 
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FIG. 2: Dependence of masses of hybrid, non-rotating com- 
pact stars on their central density for the models A, Al, and 
B. The dashed lines show the same for the associated nuclear 
equations of states. 



equations of state are rather hard. The hybrid configura- 
tions branch off from the nuclear configurations when the 
central density reaches that of the deconfinement phase 
transition. The jump in the density at constant pressure 
causes a plateau of marginal stability beyond the point 
where the hybrid stars bifurcate. This is followed by an 
unstable branch {dM/dpc < 0). Most importantly, the 
stability is regained at larger central densities: a stable 
branch of hybrid stars emerges in the range of central 
densities 1.3 < pc < 2.5 x 10^^ g cm^''. The models 
A and B feature the same high-density quark matter, 
whereas the models A and Al the same nuclear equa- 
tion of state. It is seen that the effect of having different 
nuclear equations of state (the models A and B) at inter- 
mediate densities is substantial (at densities below 10^'^ 
g cm^ all models are matched to the same equation of 
state). At the same time, the small shift 6p by which the 
models A and Al differ does not influence the masses of 
stable hybrid stars, although it is necessary for match- 
ing of nuclear and quark EOS in the models A and B. 
It is evident that there will exist purely nuclear and hy- 
brid conflgurations with different central densities but 
the same masses. This is reminiscent of the situation en- 
countered in non-superconducting hybrid stars [s^l; the 
second branch of hybrid stars was called twins, since for 
each hybrid star there always exists a counterpart with 
the same mass composed entirely of nuclear matter. 

The hybrid conflgurations are more compact than 
their nuclear counterparts, i. e. they have smaller radii 
(Fig. [3]). An exception are those conflgurations which 
are close to the bifurcation point. However, these belong 
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FIG. 3: Dependence of radii of hybrid, non-rotating compact 
stars on their central density for the model A, Al, and B. The 
dotted hnes show the same for the associated purely nuclear 
equations of states. 



to the metastable branch. Contrary to the case of self- 
bound quark stars, whose radii could be much smaller 
than the radii of purely nuclear stars, the differences be- 
tween the radii of hybrid and nuclear stars are insignif- 
icant (< 1 km) and cannot be used to distinguish these 
two classes by means of current astronomical observa- 
tions. 

Figure S] displays the astronomical bounds on the 
masses and radii of compact stars along with the "tracks" 
for our models on the mass-radius diagram. All bounds 
are quoted at the la level. The bound inferred from 
EXO0748-676, which combines information from red- 
shifted O and F lines, the emitting area of X-ray radia- 
tion, and the Eddington luminosity, constrain the mass 
and the radius of a compact star to lie on a straight line 
shown in Fig. |4] [s^ . Note that multiwavelength obser- 
vations of EXO0748-676 have led to an alternative in- 
terpretation of data which suggest a 1.35 Mq neutron 
star with a main-sequence companion (sgI . Ist} . Both the 
hybrid stars and their nuclear counterparts have their 
masses and radii within these bounds. Our sequences 
differ from those described in Refs. 0, [s^l- Instead of 
having a narrow region of stable hybrid stars, branch- 
ing off at the bifurcation point from the purely nuclear 
sequence, we obtain configurations which cover a broad 
range of central densities, comparable to those covered 
by purely nuclear equations of state. For sequences con- 
structed from models A and Al there is a range of masses 
and radii that correspond to the new family of stars dis- 
cussed in realtion with Fig. [21 The stars belonging to 
the new family lie to the left from the sharp kink at the 



FIG. 4: Mass-radius diagram for non-rotating configurations 
including the bounds from EXO 0748-676 and the bounds on 
the upper and lower pulsar masses. The sequences of hybrid 
configurations for model A (heavy, black online), Al {medium- 
light, red online), and B {light, blue online) are shown by solid 
lines. Models A and Al share the same nuclear (low-density) 
equation of state, while the models A and B share the same 
quark (high-density) equation of state. The dashed lines are 
the sequences of purely nuclear stars with underlying nuclear 
equations of state of model A (equivalently Al) and model B. 



point R ~ 13.5 km and M ~ 1.9Mq. The stable branch 
of this new family of stars is separated from the stable 
nuclear sequence by an instability region. The sequences 
corresponding to model B do not show such an instabil- 
ity region. Note that in some cases the quark equation of 
state can be indistinguishable from the nuclear equation 
of state (see, e. g., Ref. fioj) in which case quark mat- 
ter can nucleate at densities that are small compared to 
those where the purely nuclear and hybrid configurations 
separate. 

Some evidence for massive neutron stars with M ~ 
2Mq has been inferred from astronomical observations. 
A massive compact star may exist in LMXB 4U 1636-536 
with 2.0 ± O.IMq [111. Recent measurements on PSR 
B1516-f-02B in GC M5 gave M = 1.96 ± O.lMg f^. 
The models A and Al are consistent with these bounds. 
For the model B these bounds correspond to the stable 
configuration with the largest mass. 

For completeness, the lower bound on the neutron star 
mass 1.249±0.001 Mq is shown in Fig. [H which is in- 
ferred from the millisecond binary J0737-3039 [4l|. It is 
seen that the latter bound is relevant only to the low- 
central-density stars which are composed of nuclear mat- 
ter. It should be noted that canonical 1.4 Mq mass pul- 
sars will be purely nuclear if they are described well by 
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FIG. 5: Same as in Fig. [2] but for stars rotating at their 
Keplerian frequency. 

the models A and Al, but will contain quark matter cores 
if model B is more appropriate. 



B. Rapidly rotating configurations 



Millisecond neutron stars can rotate at frequencies 
which are close to the limiting orbital Keplerian fre- 
quency at which mass shedding from the equatorial plane 
starts. The Keplerian frequency sets an upper limit on 
the rotation frequency, since other (less certain) mech- 
anisms, such as gravitational radiation reaction insta- 
bilities, could impose lower limits for the rotation fre- 
quency. Rapidly rotating pulsars are potentially useful 
for placing bounds on the state of matter at high den- 
sities. Because the centrifugal potential counteracts the 
compressing stress exerted on matter by gravity, rotating 
configurations are more massive than their non-rotating 
counterparts. It is assumed that the stars are rotating 
uniformly: while nascent hybrid stars can form with large 
internal circulation, viscosity and magnetic stress will act 
to damp such motions, at least in the non-superfluid com- 
ponent [i^ l . The mass versus central-density dependence 
of compact stars rotating at the Keplerian frequency is 
similar to that for non-rotating stars (compare Figs. O 
and [5]) with the scales for mass shifted to larger val- 
ues [4^. The branch corresponding to hybrid stars is 
flattened for fast rotating configurations compared to the 
non-rotating ones, which implies that small variations in 



FIG. 6: Dependence of Keplerian (maximal) frequency of ro- 
tation for hybrid stars on their mass (in solar mass units). 
The conventions are the same as in Fig. (2] 



the stellar mass can drive the star unstable. The increase 
in the maximum mass for stable hybrid configurations 
(in solar mass units) is 2.052 2.462 for the model A, 
2.017 2.428 and 2.4174 for the model Al (there are 
two maxima) and 1.981 2.35 for the model B. 

Fast spinning radiopulsars can potentially limit the 
mass-radius relation of hybrid stars. The requirement 
that the maximum spin frequency cannot exceed the 
Keplerian frequency at the neutron star surface trans- 
lates into the bound fmax < 10A5{M/Mq)'^/'^R~q^^, 
where R is the radius of non-rotating star in units of 
10 km [H,!!^. The frequency of the fastest ra diop ulsar 
observed to date, PSR J1748-2446ad, is 716 Hz The 
dependence of the orbital Keplerian frequencies of hybrid 
stars on their maximum mass is shown in Fig. [HI The con- 
figurations with masses below 1 Mq are excluded by the 
observed frequency 716 Hz. However, hybrid stars are 
much more massive - they bifurcate from the sequence of 
purely nuclear configurations for M > 2Mq. Therefore 
this particular observation is not useful in limiting the 
properties of CCS hybrid stars. 

Finally note that at the point of bifurcation the Keple- 
rian frequency of configurations jumps to higher values, 
which is the consequence of the fact that hybrid stars, 
having the same mass as their nuclear counterparts, are 
more compact and thus can support larger rotation rates. 

The discovery of the double-pulsar system PSR J0737- 
3039 [111, [U offers a unique opportunity to place further 
bounds on the gross parameters of compact stars by a 
measurement of the moment of inertia of star A, since the 
spin frequencies and the masses of both pulsars are ac- 
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FIG. 7: Dependence of moment of inertia of the hybrid con- 
figurations on the mass of configuration (in solar mass units) 
rotating at the limiting frequency. The conventions are the 
same as in Fig. [2l 



curately measured. Timing measurements over a period 
of years could provide information on spin-orbit coupling 
which could be revealed through an extra advancement 
of the periastron of the orbit above the standard post- 
Newtonian advance or in the precession of the orbital 
plane about the direction of the total angular momen- 
tum [4^. The dependence of the moment of inertia, /, of 
configurations on their mass (in the case of rotation at the 
Keplerian frequency) is shown in Fig. [T) Since the mo- 
ment of inertia is independent of the rotation frequency, 
we have chosen to extract it for configurations rotating at 
the limiting frequency; note that the masses of configura- 
tions should be rescaled appropriately, if one is interested 
in the I(M) function for slowly rotating configurations 
(see ref. [49[ for a computation of the moment of inertia 
in the slow-rotation approximation and a discussion of 
the moment of inertia in the context of PSR J0737-3039). 
It is seen that while an accurate measurement of the mo- 
ment of inertia of pulsar A can discriminate between the 
two nuclear equations of state, it will not be useful for 
accessing the properties of hybrid stars since the mea- 
sured masses of pulsars are too low: M/Mq = 1.337 
for pulsar A and M/Mq = 1.250 for pulsar B; as ap- 
parent from Fig. [7] observations of heavier objects are 
needed to obtain useful bounds on the moment of inertia 
of hybrid configurations. The differences ~ 20% in the 
moments of inertia of purely nuclear and hybrid config- 
urations of the same mass (in the case of the models A 
and Al) are within the accuracy that can be achieved 
in measurements similar to PSR J0737-3039. Finally we 
note that the moment of inertia of the CCS phase is re- 



quired to compute the strain amplitude of gravitational 
wave emission from an isolated spinning CCS hybrid star. 
Upper limits on the strain amplitude have been placed 
by the LIGO and GEO600 detectors HH. The key 
unknown quantity is the equatorial ellipticity of a hy- 
brid star e = [I^x — Iyy)/Izz, where la are the compo- 
nents of the moment of inertia. The latter quantity has 
been estimated for CCS stars at the qualitative level in 
Ref. [53 | assuming uniform-density incompressible mod- 
els and maximally strained core [53| . Our models of CCS 
stars can be utilized to obtain quantitative limits on the 
strain of gravitational wave emission that can be emitted 
by CCS hybrid stars. 



IV. CONCLUSIONS 

We constructed a set of equations of states of hybrid 
stars which are composed of color-superconducting quark 
matter at high, and purely nuclear matter at low densi- 
ties. High-density quark matter is described in terms of 
the semi-microscopic NJL model which includes pair cor- 
relations that lead to the three-flavor LOFF phase as the 
ground state of QCD at moderate densities. The low- 
density nuclear phase is described in terms of hard rela- 
tivistic equations of state based on the Dirac-Bruckner- 
Hartree-Fock theory. The sequences of stellar configura- 
tions constructed on the basis of these equations of state 
reveal a new branch of stable hybrid configurations featur- 
ing CCS matter. These are analogous to the twin config- 
urations obtained for non-superconducting quark matter 
in Refs. 34] and are in contrast to the previously reported 
NJL-model based sequences of non-rotating three-flavor 
color-superconducting hybrid stars, which were either en- 
tirely unstable or had a narrow range of stable configu- 
rations adjacent to the point where they bifurcate from 
the nuclear ones. 

A further search in the parameter space is required 
in order to judge how generic our constructions are; a 
quite general conclusion of our analysis is that the nuclear 
equations of states need be hard to enable the matching 
to the quark equations of state. Furthermore, matching 
is only achieved if the ratio of the couplings in the diquark 
and quark-anti-quark channels is large (77 = 1) compared 
to the value obtained through the Fierz transformation 
(77 = 0.75). The lack of confinement in the NJL model 
leaves the low-density behavior of the quark matter equa- 
tion of state somewhat uncertain. Small variations in the 
bag constant help in achieving a phase transition, but do 
not affect largely the global properties of hybrid stars 
(cf. our models A and Al). The maximum masses of 
our stable hybrid configurations are consistent with the 
presently available astronomical bounds on masses and 
the mass-radius relationship of compact stars. 

We provided a first discussion of rapidly rotating con- 
figurations of CCS hybrid stars and mapped out their 
masses and orbital Keplerian frequencies at which they 
are destabilized by mass-shedding from the equator. 
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Rapidly rotating configurations support larger masses 
(as expected). The mass-central density relation for this 
class of objects is "flatter" than for their non-rotating 
counterparts, i. e., smaller variations of mass can drive 
the system unstable. The Keplerian frequencies of mas- 
sive members of the sequences, in particular hybrid con- 
figurations, are by a factor of two larger than the rotation 
frequency 716 Hz of the fastest millisecond pulsar known 
to date. We find that the hybrid configurations, being 
more compact, can rotate at substantially larger rates 
than their nuclear counterparts of the same mass. The 
moment of inertia of a hybrid star is smaller than that 
of a nuclear star of the same mass. At the point of bi- 
furcation it drops abruptly giving rise to a sequence of 
stable hybrid stars with almost the same masses but dif- 
ferent moments of inertia; this implies that simultaneous 
mass and moment-of-inertia measurements are useful in 



distinguishing hybrid configurations from their nuclear 
counterparts. 
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